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AN ANALYTICAL STUDY OF THE EFFECT OF AIRPLANE WAKE ON THE LATERAL DISPERSION 
— OF AERIAL SPRAYS ' 

By WiLMEB H. Reed III 


SUMMARY 

Calculations are made to determine the trajectories of liquid 
droplets introduced into the air disturbances generated by an 
airplane engaged in aerial spraying. The effects of such factors 
as the position at which droplets are ejected into the disturbances, 
airplane lift coefficient, and altitude are investigated. The 
distribution of deposit on the ground is computed for several 
droplet-size spectra, variations in the rate at which mass is 
ejected along the span, and lateral flight-path spacings. Con- 
sideration is then given to the problem of adjusting these factors 
with the aim of improving the uniformity and increasing the 
effective width of the deposit. 

The results indicate that the lateral dispersion of droplets is 
increased when the spanwise position at which particles are 
ejected is moved toward the wing tip. Greater dispersion also 
results when the airplane lift coefficient or altitude is increased. 

With the spray discharged continuously along the span at a 
constant rate, the deposit has a maximum concentration at the 
plane of symmetry and rapidly diminishing concentration 
beyond the wing tip. In such cases, it was found that the uni- 
formity and effective width of the swath could be improved by 
increasing the mass effiux rate with spanwise distance from the 
plane of symmetry. When this is done, the lateral distance 
between adjacent flight paths required for a given degree of 
uniformity shows an increase amounting to from O .4 to 1.0 
semi-span. Of the two droplet-size spectra considered, the one 
having the smaller mean diameter {200 microns) gives a more 
uniform deposit for a given flight-path spacing and the degree 
of uniformity is less sensitive to changes in the spacing of 
adjacent passes of the airplane. 

INTRODUCTION 

Although airplanes have been used to disperse liquid 
sprays and dust for insect control and agricultural purposes 
for many years, only recently has an appreciable amoimt of 
study been devoted to the problem of improving the efficiency 
and effectiveness of the operation. One of the problems 
recognized as being important is that of obtaining a suitable 
distribution of the deposit. Experimental studies made by 
the Forest Insect Laboratory of BeltsviUe, Maryland (ref. 1) 
have indicated that, in general, aerial spray deposits exhibit 
excessive peaks of concentration beneath the airplane with 
diminishing intensity toward the outer edge of the swath. 


Although it is possible to improve these measured deposit 
characteristics by a tedious trial-and-error process, it is 
believed that a theoretical study of the important physical 
factors involved would provide a better understanding of the 
mechanism which determines the lateral dispersion of parti- 
cles ejected into the flow field behind an airplane. Such a 
study may also suggest possible methods of improving the 
deposit characteristics which might otherwise be overlooked. 

In the present report the paths of liquid spray droplets 
issued into the flow field behind an airplane are calculated. 
Considered in these computations are the effects of such 
factors as the droplet-size spectrum, spanwise location of 
nozzles, height of airplane above the ground, and aero- 
dynamic flow field in the wake. Distributions of deposit 
are determined for several variations of these parameters, 
with a view to improving the uniformity and the effective 
width of the swath. The airplane used as an illustration 
in the analysis is the agiicultural prototype designated 
the Ag-1. 

SYMBOLS 

A wing aspect ratio 

b wing semispan, ft 

C concentration of deposit, Ib/unit area 
Cd drag coefficient of droplet 

Ci lift coefficient of airplane 

D drag force on droplet, lb 

g acceleration due to gravity, ft/sec^ 

m mass of droplet, slugs 


Vi—yo 

Q mass ejected per unit length along flight path, Ib/unit 
length 

q velocity induced by isolated rectilinear vortex, ft/sec 

R Reynolds number of droplet 

Ru Reynolds number of droplet moving through still air 
at velocity U 

r radial distance from vortex center, ft 

t time, sec 

U flight velocity, ft/sec 

V percent of total volume 

»o,Wa velocity components of air in y- and s-direction, re- 

spectively, ft/sec 


J Supersedes NACA TN 3032, “An Analytical Study of the Effect of Airplane Wake on the Lateral Dispersion of Aerial Sprays” by Wilmer H. Reed III, 1053. 
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Vd,Wa 


s' 


OL 


r 

5 




M 

V 

Pa 

Pd 

a 

T 


velocity components of droplet in y- and s-direction, 
respectively, ft/sec 

longitudinal, lateral, and vertical coordinates of drop- 
let, respectively, semispans 
vertical coordinate of droplet with origin shifted to 
the bound vortex, semispans 


dimensionless parameter. 


ffCoRt 

~2^ 


circulation, ft^/sec 

diameterof droplet, microns (1 micron=3.28X10“®ft) 
mean diameter of droplets, microns 
dispersion of droplet diameters, microns 
lateral and vertical coordinates of trailing vortex 
cores, semispans 

absolute viscosity of air, slugs/ft-sec 
Idnematic viscosity of air, ft^sec 
density of air, slugs/cu ft 
density of droplets, slugs/cu It 

dimensionless parameter, 18 | — 

0 Pi 

dimensionless time, C/i/6 


Subscripts: 

0 initial value 

g intersection of droplet path with ground plane 

y,z lateral and vertical directions, respectively 

t terminal velocity of droplet 

c correction due to bound vortex 


EQUATIONS OF MOTION 

The motion of spray droplets in a space-fixed plane located 
behind the airplane perpendicular to the line of flight will 
be considered. A rectangular-coordinate system is used 
where vertical z distances are measured in airplane semi- 
spans from the ground plane (positive upward) and lateral y 
distances in semispans from the line of flight. If Vi and Wi 
are, respectively, the lateral and vertical components of 
droplet velocity and Va and Wa are the correspondiug com- 
ponents of air velocity in the plane,, the equations of motion 
which establish dynamic equilibrium of the forces acting on 
a droplet are 


dvd 7 -> 

~-g=D, 

(1) 

dWi „ 

(2) 


The resultant aerodynamic drag force on the droplet is 

and the lateral component of drag which appears in 
equation (1) may be expressed 


■\KVa—Vif+{Wa — Wif 


ANALYSIS 

The paths taken by spray droplets issued from an airplane 
are influenced by the flow field induced in the region of air 
through which the airplane has flown. When drag forces 
acting on the droplets are predominant, as in the case of 
small particles, the trajectories follow the streamlines of 
flow more closely than do the trajectories of larger particles 
for which inertia forces are predominant. In aerial sprays, 
the smallest particles may remain in suspension in the air 
indefinitely, whereas the largest particles fall in paths prac- 
tically uninfluenced by the surrounding air flow. In this 
section the method of determining trajectories of aerial 
spray droplets and the resulting distribution of deposit is 
presented. 

Calculations of a similar nature have been performed by 
various investigators of aircraft icing problems. In these 
studies the trajectories of small water droplets moving in 
an airstream past various bodies such as airfoils and cylinders 
have been calculated to determine area and rate of droplet 
impingement on the surface. Glauert, in reference 2, con- 
sidered cases wherein the air velocities and droplet diameters 
were sufficiently small so that the drag could be 

assumed proportional to the velocity of the droplet relative 
to the air (Stokes’ law of resistance). Later, Langmuir and 
Blodgett (ref. 3) extended this work and considered the 
variation of drag beyond the limits of velocity and diameter 
in which Stokes’ law is valid. It was found, after a considera- 
tion of the range of droplet Reynolds number encompassed 
in aerial spraying, that the deviation of drag from Stokes’ 
law should be included in the present calculations. 


According to Stokes’ law (see ref. 4, p. 598) the drag of 
spheres at very low Reynolds numbers (i?<^l) is 

^ 

i? Stokes = 67T/i — 

Now, since 

^ _ PaS ^|{Va — Vdy+ {Wg — W^^ 

M 


the ratio of the actual drag to the drag based on Stokes’ 
law becomes 

D CqR 

17 stokes 24 

Since 



and 


t) Pa^U 

Ku— 

p 



the nondimen sional equation of motion in the ^/-direction, 
after appropriate substitutions into equation (1), becomes 

diVdjU) a- CJR ( Va Vi\ 

dr Ra 24 \U U) ^ ’ 

where 

tr=is|- 

0 Pd 

In a similar manner the equation of motion in the z-direction 
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(eq. (2)) becomes 

d{wa[U) a CdR /Wa wA bg 

dr Rv 24 \U U) ^ ^ 

Except for the inclusion of the term due to gravity, 
bglU^, which appears m equation (4), equations (3) and (4) 
are of the sam e form as the equations of motion given in 
reference ST^^Since the left-hand sides of equations (3) and 
(4) represent the accelerations of the droplet in the lateral 
and vertical directions, a double integration of these equa- 
tions with prescribed initial values yo, So, ivalU)o, and (w<r/f7)o 
defines the subsequent positions of the droplet. 

In carrying out the integration of equations (3) and (4) 
it is necessary to evaluate R and CnRj24: as the particle 
progresses along its trajectory. The Reynolds number of the 
droplet is computed from the relation 

/i?Y /Va VaV /Wa WaV 

^yu-u) +yu-u) 

and from reference 3 the term CdRI24: is tabulated as a func- 
tion of R in table I. 

AIR VELOCITIES 

The air velocities Va and Wa as functions of droplet position 
and time are needed to solve equations (3) to (5). Air 
disturbances which are not related to the airplane, such as 
atmospheric turbulence and wind, are not dealt with in the 
present analysis. In the wake of an airplane the major 
disturbances are associated with the propeller wash and the 
trailing vortices of the wing. Of these, the effect of propeller 
wash is largely local! ed so the vicinity behind the fuselage 
and its air disturbance velocities are approximately in the 
direction of the flight path. Because the inclusion of these 
and other air velocities in the direction of flight would 
greatly complicate the analysis and since such disturbances 
are believed to have a small effect on the lateral displacement 
of droplets, they have been neglected. The wing, on the 
other hand, sheds a sheet of vortices which soon roll up into 
two discrete vortex cores that are initially spaced somewhat 
less than a wing span apart. These so-called trailing vortices 
and the velocities induced by them in the i/ 2 -plane persist for 
a considerable time after the airplane’s passage and appear 
to be of fundamental importance in determining the lateral 
distribution of droplets which impinge on the ground. 

Idealized flow pattern. — An idealized flow pattern, void 
of viscosity, is used herein to approximate the velocities 
which are induced behind a wing. The trailing vortices 
and their images reflected from the ground plane are repre- 
sented by a system of four parallel rectilinear vortices 
which pass perpendicularly through the j/s-plane at points 
A, B, B', and A' (see fig. 1) and extend to infinity ahead of 
and behind this plane. The vortices are of equal strength 
and have the directions of rotation indicated by the figure. 
The coordinates of the vortex cores of this system change 
with time, since the coordinates of one of the vortex cores 
moves at the velocity induced at that point by the remain- 
ing three vortices. As a resuljt, the trailing vortices are 
displaced downward and outward as depicted by the dashed 
fines in figure 1 . 



Figure 1. — Paths taken by the trailing vortices in the presence of a 

ground plane. 


If the equations which express the velocities of the cores 
and the coordinates of their paths (ref. 4, pp. 223 to 224) 
are combined and integrated with respect to time, the 
positions of the cores as functions of time are obtained. If 
a rectangular spanwise loading is assumed, the strength of 
the trailing vortices is expressed as 

r^C^Ub 

A 

and these relations then become 


Cr.0, 




( 6 ) 

( 7 ) 


where and f are, respectively, lateral and vertical coordi- 
nates of the vortex core in semispans and the constants of 
integration are 


1)0 So 


^ ClT]!? 2 

p 2 


The velocity components vJU and wJU induced by the 
vortex system at a point having coordinates y and z semi- 
spans are, by the Biot-Savart theorem. 


Va Cl r f +2 f +2 

U~2tA hi+zf+iv-y? A+zf+iv+yf 


Wa Cl r v—y I vAy 

U~2TA\_AAzYA{r)-yY'^AAzYA{r,AyY 

vAy v—y ~| 

A-^YA(vAyY A-^YAA-vYA 


( 8 ) 


( 9 ) 
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The streamlines indicating the direction of air velocity 
vectors at a given time are shown in figure 2. The effect 
of foliage and other obstructions to air flow near the ground 
is, of course, neglected by the idealized flow pattern. 

Because the induced velocities given by equations (8) 
and (9) become infinite at the position of the line vortices, 
it is necessary to consider the minimum radial distance 
from these points at which these equations remain appli- 
cable to physical conditions where infinite velocities are 
not possible. In reference 5 an approximate expression is 
derived for the radius of the rolled-up vortex cores behind 
a wing. It is shown by simple energy considerations that 
the core radius for a wing having an elliptical spanwise 
loading is 0.155 semispan. In order to avoid the region of 
the vortex core, the spanwise starting positions of the 
droplets should be at least this distance from the center of 
the core. 

Velocities in a real fluid. — Since the determination of 
velocities in the assumed flow field is based on an ideal fluid, 
it is of interest to consider briefly the effects of viscosity on 
the flow velocities in a real fluid, or, more specifically, the 
dissipation of vorticity with time. 

To illustrate the effect of the decay of vorticity on the 
velocities induced in a real fluid, consider an isolated recti- 
linear vortex with an initial strength of To. At a subse- 
quent time t, the circulation in a circle of radius r about the 
vortex (see ref. 4, p. 592) is 

r = ro ( l - e -^'^^^0 (10) 

Hence the tangential velocity on the circle becomes 

( 11 ) 


Equation (12) is plotted against r for various values of time 
in figure 3. The figure shows the extremely slow rate at 
which induced velocities a few feet from the vortex line are 
diminished as a result of viscosity. An additional viscous 
effect not accounted for by the idealized flow pattern is 
evidenced in the form of a boundary layer of reduced air 
velocities adjacent to the ground. 

TRAJECTORY COMPUTATIONS 

The sample airplane for which the trajectory computa- 
tions were made is the Ag-1, a single-engine, low-wing, 
agricultural airplane of moderate si e. The liquid droplets 
are assumed to have the density of kerosene, a frequently 
used solvent for insecticides. The physical characteristics 
pertinent to these computations are given in table II. 

A step-by-step integration of equations (1) and (2) was 
accomplished on the Bell Telephone Laboratories X-66744 
relay computer at the Langley Laboratory. It was found 
convenient to tabulate the variation of CdB ! 24i with R for 
several values of R and use a linear interpolation between 
the tabulated values. The coordinates of the trailing vortex 
cores were determined in a similar manner by tabulating 
T) and f at suitably chosen values of r. The air velocities 
were computed at each point along the trajectory by eval- 
uating equations (8) and (9) directly. 

A series of trajectories was computed for several flight 
conditions with selected values of the parameters o- and Ru 
and with assigned initial conditions j/o, ^o, (®tf/L)o, and 
{WdlU)o. In order to reduce the number of computed 
points along a trajectory and yet retain a reasonable degree 
of accuracy, the time increments during an integration were 
varied; the criterion used to select these time intervals was 


and therefore 


Q_viscous 

Q_ideal 




( 12 ) 


z 



Time, min 



Figure 3.— Effect of viscosity oil the reduction of velocity induced by 
an isolated rectilinear vortex as a function of time and distance from 
the core. 


Figure 2. — Streamlines induced by the trailing vortices in the presence 
of a ground plane. 
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the change in the droplet acceleration evaluated at two 
adjacent points. The integration was stopped when the 
particle intersected the ground plane or when its trajectory 
looped about the vortex core. 

The vortex cores are assumed initially to be at the height 
of the wing and separated laterally a distance of 1 span. 
Particles are given starting positions at various points along 
a line connecting the trailing vortices, and the initial drop- 
let velocity is taken to be the steady-state velocity of the 
droplet if it were falling under the influence of gravity in a 
uniform airstream having velocity components equal to the 
induced velocities of the idealized flow pattern at the point 
of release. This assumption was made, rather than the 
assumption that the particles were released from rest at 
the trailing edge, in order to compensate somewhat for the 
error caused by neglecting the bound vortex in regions close 
to the wing (see appendix A). Actually, in the neighbor- 
hood of the wing, air velocities are induced by a bound 
vortex at the wing and, in addition, a continuous sheet of 
vorticity extending downstream from the bound vortex. 
In the present case, where a simplified uniform spanwise 
loading is considered, the flow in this region may be repre- 
sented by a single horseshoe-vortex arrangement with a 
bound vortex at the quarter-chord point connected with 
trailing vortices of equal strength which terminate at the 
wing tip and extend downstream to infinity. By neglecting 
the bound vortex, the initial downwash is underestimated, 
although this error in downwash is to some extent reduced 
by the forward-extending vortices assumed in the idealized 
flow pattern. 

At a distance of 2.0 semispans behind the bound vortex 
the velocities induced by the idealized flow pattern were 
found to be within 6 percent of the velocities induced by 
the horseshoe-vortex system. The approximate error in 
the vertical position of a droplet 2.0 semispans downstream 
of the wing, caused by an underestimation of the down- 
wash ahead of this point together with the compensating 
assumption of a finite downward ejection velocity, is eval- 
uated in appendix A. 

A summary of the results of droplet-trajectory computa- 
tions is given in table III. Here, the quantities yo, yg, 
u>a^, and tg are given for the various droplet sizes and airplane 
flight conditions considered in the analysis. The calculated 
droplet trajectories for these cases are plotted in figure 4. 

Figure 4 indicates that as the droplet diameter is reduced 
at a given nozzle location the influence, of disturbance veloci- 
ties becomes more pronounced, until finally a critical condi- 
tion is reached wherein the path of the droplet encircles the 
vortex core. It appears from the limited number of trajec- 
tory computations available that once the particle encircles 
a vortex core it continues to do so for a considerable time, 
during which it is carried laterally outward with the core. 
Inasmuch as the labor involved in calculating the complete 
paths of such particles would be excessive and, owing to 
extensive scattering, the corresponding concentration of 
deposit would probably be very small, droplets which 
initially loop about the vortex eenters are assumed to. remain 
suspended in the flow field indefinitely. The percentage of 


mass represented by these droplets (mass not recovered in 
the deposit) is evaluated later in the report. 

An additional effect not accounted for in the analysis is 
the evaporation of droplets during descent. It might be 
noted that, since, in general, droplets which encircle the 
vortex cores remain aloft for the greatest time, evaporation 
losses are probably most predominant in this pai-t of the 
spray. Evaporation effects, therefore, tend to compensate 
for the previous assumption which involved the omission 
of part of the deposit. 

DISTRIBUTION OF DEPOSIT 

From the foregoing trajectory computation it is possible 
to determine the distribution of deposit for various assumed 
nozzle characteristics and spanwise rates of efflux. Consider 
first the case where a spray, composed of particles of various 
sizes, is introduced into the flow field at point (i/o, So) and at 
the rate of Q pounds per unit length along the flight path. 
For each droplet size there is a corresponding position of 
impingement on the ground yg; therefore, droplets of size d 
and 5-fA5 fall at yg and yg-\-Ayg, respectively. The rate at 
which mass is deposited in element Ayg is Q Av where Av is 
the percentage of Q which is represented by droplets having 
diameters between 5 and 5-fA5. The concentration of 
deposit at position yg is then 


C=Q 


Av 


which, in the limit, can be written 


^ ^ dv d8 


(13) 


(14) 


Equation (14) applies for a single nozzle, and by super- 
position the concentration of deposit can be determined for 
any number of nozzles. For a- large number of nozzles 
closely spaced along the wing, however, it is more convenient 
to consider mass as being ejected, cohtinuously along the 
span rather than at a finite number of discrete points. 

The rate of change of concentration at yg with respect to 
Q is, by equation (14), 


dC dv dS 

dQ d8 dyg 


(15) 


but, since is a continuous function of j/o. 




Therefore, the total-concentration at yg due to a contmuous 
mass flow along the span becomes 


r‘'o dQ & d^ 

Jo dyodddyg 


dyo 


(16) 


In order to compute the concentration of deposit by means 
of equation (14) or (16), suitable expressions for the terms 
dv/d& and dbjdyg must be determined. The relation dvjdb is 
defined by the nozzle or spray characteristics, and dbjdyg is 
associated with the trajectories of droplets. 









y, semispans 


(b) Ct=2.2; 2 o= 0.5 semispan. 


Figure 4. — Continued. 


325362—55 2 




semlspans _ semispans 


8 


BBPORT 1196 — NATIONAL ADVISORT COMMITTEE FOR AERONAUTICS 



y, semispons 


(c) Ct= 1.2; Sb= 1.0 semispan. 
Figure 4. — Concluded. 
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NOZZLE AND SPRAY CHARACTERISTICS 

Reference 6 shows that the breakup of liquid jets, ejected 
into a high-velocity airstream, is dependent on the Reynolds 
number of the mean droplet size and the relative viscosities 
of the liquid particles and the surrounding air. The nozzle 
bore was found to have no effect on the mean droplet size 
except at very low ejection velocities, in which case the 
droplets were more uniform and had diameters approxi- 
mately twice that of the nozzle. The mean droplet diameter, 
defined as 8„, is considered herein to he such that half of the 
total mass contained in the spray is represented by droplets 
smaller than 5m and the other half, by droplets larger 
than Sm- 

An empirical expression for the mean diameter, which is 
applicable for all conditions lilcely to be encountered in aerial 
spraying, is given in reference 6. In the notation of the 
present report the expression is 



Figure 5.— -Sample droplet-size spectra used in analysis. 


where 5^ is the mean droplet diameter in cm, V is the relative 
velocity in cm/sec between the liquid jet and the air, and Vd 
is the kinematic viscosity of the liquid in cmVsec. 

Inasmuch as the droplet-size distribution curves measured 
in reference 6 were nearly symmetrical about Sm, the so- 
called Gaussian or normal distribution gives a very good 
representation of the spectra. The equation of the normal 
distribution (see ref. 7) then provides the following mathe- 
matical expression for the term dv/d8 appearing in equations 
(14) and (16): 


r o.675(;-;„) n 

dv^ Q.Qld L « J 
dS e-y]2ir 


(18) 


The term e in equation (18) is known as the “probable 
error” and in the present problem represents the dispersion 
or spread of droplet sizes about 6^/ that is, e is such that 50 
percent of the volume is composed of droplets having diam- 
eters between 6^—6 and Sm+e- The magnitude of e is then 
a measure of the uniformity of droplet sizes in a liquid spray 
having a mean diameter Sm- 

Two sample spectra of droplet sizes are considered in the 
present analysis; these have mean diameters of 200 and 300 
microns and a dispersion e of 50 microns. The percentage of 
the total mass in terms of droplet diameter may be de- 
termined by integrating equation (18)., a plot of which is 
given in figure 5 for the two droplet spectra under considera- 
tion. (See tabulated values of the probability integral in 
ref. 8.) 

MATHEMATICAL REPRESENTATION OF TRAJECTORY DATA 

An empirical relation which expresses the droplet diameter 
as a function of yo and yg for a given flight condition is 


derived in appendix B by applying the method of least 
squares to the trajectory data presented in table III. An 
equation which was found to give a satisfactory approxima- 
tion of the computed data is 

S = ki-\-k2P-{-k^P^ (19) 


where 


Vs—yo 

and the coefficients fc], k 2 , and kg are functions of yo- It is 
seen that as yg becomes large the droplet diameter ap- 
proaches the limiting value k^. The mass represented by 
droplets of diameter ki and smaller corresponds to the mass 
previously omitted on the assumption that small droplets 
whose paths encircle the vortex centers are entrained by the 
air. 

On this basis the percentage of mass not recovered is 
illustrated in figure 6 for the two sample spectra with 
Ci,= 1.2 and 2 q= 0.5 semispan. 



Figure 6. — Percentage of mass not recovered in deposit as a function 
of spanwisB starting position for Ci,= 1.2 and Zo— 0.5 semispan. 
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RESULTS AND DISCUSSION 

DISTRIBUTION OF DEPOSIT WITH DISCRETE NOZZLES 

When liquid is ejected behind the wing from a discrete 
point, such as a nozzle, the resulting concentration of deposit 
on the ground may be determined by evaluating equation 
(14) with the aid of equations (18) and (19) at various yg 
positions on the ground. The distribution of deposit so 
computed is presented in figures 7 to 9 for a unit mass efflux 
rate. These figures are intended to illustrate the effect of 
spanwise location of the nozzle, mean diameter of the drop- 
lets, airplane hft coefficient, and airplane height. 

The effect of the mean diameter of the droplets on the 
shape of the deposit curves is shown in figure 7 for the case 
of (7i=1.2 and Zo=0.5 semispan. Here, the deposit curves 
are characterized by high peaks of concentration when the 
nozzles are located near the plane of symmetry (i/o=0.25 
semispan), whereas for nozzle locations in the vicinity of the 
wing tip ( 2 / 0 = 0.75 semispan) the concentration of deposit is 
less intense and the distribution is more uniform. This 
effect would be expected since the magnitude and curvature 
of the air-flow velocities are increased as the wing tip is 
approached, and, as a consequence, droplets of different 
diameters are dispersed more in this region than at points 
farther inboard. 

Curves for various values of mean diameter of the droplets 
indicate that as is reduced the deposit is distributed over a 
larger area, the center of which is displaced slightly farther 
outboard. 

The effect of airplane lift coefficient on the deposit curves is 
illustrated in figure 8, where the altitude is again taken to be 
0.5 semispan and 5^=300 microns. At (7x,=2.2, the lift 
coefficient corresponding to the flaps-down condition, the 
deposit curve is of similar shape but displaced farther out- 
board than for the ffaps-up case (<7z,= 1.2). It is of interest 
to note that the velocity of impingement for the flaps-down 
case (see table III) is generally slightly higher than for the 
corresponding flaps-up case. The higher droplet velocities 
and also the stronger air-disturbance velocities which accom- 
pany the use of flaps may provide greater penetration of the 
crop foliage by the sprayed liquid. 



Figube 7. — Effect of the mean droplet diameter on distribution of 
deposit for discrete nozzles. Ci,= 1.2; zo=0.5 semispan. 


A comparison of the deposit curves calculated for flight 
altitudes of 0.5 and 1.0 semispan is shown in figure 9, where 
f7i=1.2 and 5^=300 microns. The figure indicates that the 
uniformity and lateral displacement of the deposit is in- 
creased with airplane altitude. 

DISTRIBUTION OF DEPOSIT WITH CONTINUOUS SPANWISE MASS FLOW 

In practice, aerial sprays are generally issued from many 
nozzles along the span. In order to calculate the total 
concentration of deposit when there is a multitude of closely 
spaced nozzles, it is convenient to use equation (16), which 
is the integral form of equation (14). Variations of dQjdy^ 
were assumed and the resulting concentrations of deposit 
were determined by integrating equation (16) numerically 
with the use of Simpson’s rule. In performing these com- 
putations the available trajectory data (table III) were 
extrapolated by means of equation (19) so as to encompass a 
range of spanwise starting positions from 2 / 0 = 0 to 2 / 0 = 0.8 
semispan. In the results which follow, Cx,= 1.2, 2 q= 0.5 
semispan, and the mass efflux rate is assumed to be zero 
beyond 2 / 0 = 0.8 semispan. 



Figure 8. — Effect of airplane lift coefficient on distribution of deposit 
for discrete nozzles. Zo = 0.6 semispan; 6„=300 microns. 



Figure 9. — Effect of airplane height on distribution of deposit for 
discrete nozzles. Cl= 1.2; a„=300 microns. 
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In figure 10 are presented deposit curves evaluated for a 
uniform spanwise efflux rate (solid lines in the figure) and 
also deposit curves in which dQldy^ was modified in an attempt 
to improve the uniformity of deposit (dashed lines). For 
the case of a uniform spanwise mass efflux rate the deposit 
curves for both of the droplet-size spectra considered show 
maximum concentration at the plane of symmetry. The de- 
posit curve' for 5^=300 microns (fig. 10 (b)) is the less 
uniform of the two ; it has a higher concentration at the plane 
of symmetry but at yg=2.0 semispans it has about one-half 
the concentration of the deposit curve for 3^=200 microns 
(fig. 10 (c)). The deposit is, of course, reduced to zero in all 
cases as y* becomes large. 

For the purpose of determining a modified spanwise efflux 
rate which would improve the uniformity of deposit, assume 
that all of the droplets issued from the wing were of equal 
size (S=Constant) . Each of the droplets ejected at a point 
yo would then fall at a single point y^ on the ground in ac- 
cordance with the trajectory characteristics of the droplet 
size in question. The concentration of deposit in this case 
would be 


Q^dQ dyo 
dyo dyg 


( 20 ) 


Therefore, the variation of dQ/dyo with required to provide 
a uniform deposit of a constant-droplet-size spray becomes 

Constant 
dyo~ dyo/dyg 

where dyoldyg is expressed as a function of yo alone by the 
simultaneous solution of equation (19), with 6 equal to a 
constant, and the relation 

dyo_ bd/byg 

dyg dS/dyo ^ 

Although this procedure for determining the optimum span- 
wise variation of mass flow rate is true only for sprays of 
constant droplet size, the uniformity of deposits having a 
variable droplet size should be improved by varying the mass 
flow rate in accordance with equation (21), taking the droplet 
size to be the mean diameter of the droplets. 

The modified efflux rate I in figure 10(a) is equation (21) 
evaluated with 5=300 microns. The constant in equation 
(21) was selected to give the same total mass flow per unit 
length along the flight path as was assumed for the case of 
a constant spanwise flow rate (Q=0.8 pound per unit length). 
The resulting deposit is shown by the appropriate coded curve 
in figure 10(b). Here, the peak deposit which occurs at 
y,= l.l semispans indicates that too large a portion of the 
mass is transported laterally outward. One method of over- 
coming this difficulty would be to evaluate equation (21) for 
a droplet diameter somewhat larger than 5„. 

The mass flow rate evaluated to give a uniform deposit 
with a constant-diameter spray of 350 microns is denoted 
by II in figure 10(a). The corresponding deposit for the 
case of a spectrum of 300-micron mean diameter is shown 
in figure 10(b). Here, the distribution of deposit is seen to 
be approximately uniform up to yg=1.2 semispans, and be- 
yond this point it breaks off abruptly as do the other deposit 
curves in figure 10(b); however, the lateral extent of the 


deposits of variable mass flow rate is somewhat greater than 
the deposit calculated for a uniform mass flow rate. 

Since a more favorable deposit results when the diameter 
used to evaluate equation (21) is somewhat larger than the 
mean diameter of the droplet spectrum, efflux rate I was 
used to modify the deposit curve for 5„=200 microns. The 
resulting deposit is shown in figure 10(c). The concentra- 
tion of the modified deposit for 6^=200 microns is less in- 
tense but is uniform over a larger region than the modified 
deposit with 6„=300 microns. Moreover, there is no abrupt 
reduction in the intensity of concentration beyond y*=1.2 
semispans as was the case for 6^=300 microns. 

In addition to varying the mass flow rate along the wing, 
there remains the possibility of regulating the mean diameter 
of the droplet spectrum at various nozzle locations along the 
wing. As indicated in equation (17), this could be accom- 
plished through control of the relative jet velocity, which is 
dependent on the direction of the jet relative to the smround- 
ing airstream and the pressure under which liquid is ejected. 

To investigate one such case, it was assumed that the mean 
diameter increased linearly from 200 microns at the plane of 




(a) Mass efflux rate variations. 

(b) Deposit curves; 300 microns. 

(c) Deposit curves; 5„=200 microns. 

Figube 10. — Effect of spanwise variations of a continuous mass efflux 
rate on distribution of deposit. Ci, = 1.2; Zo=0.5 semispan. 
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symmetry to 300 microns at yo=0.8 semispan and that the 
spanwise efflux rate was held constant. In this example, 
the spectrum of smallest mean diameter (5=200 microns) 
was issued from the plane of symmetry in order to reduce the 
peak deposit which occurred in this region when sprays of 
similar spectra were issued at a uniform rate across the span. 
The distribution of deposit for this case is plotted in figure 11. 
It is seen that this deposit curve has characteristics similar 
to the modified deposit of figure 10(b), where the uniformity 
of deposit was improved by increasing the mass fiow rate with 
distance from the midspan. 

EFFECT OF FlIGHT-PATH SPACING ON DISTRIBUTION OF DEPOSIT 

Thus far, an attempt has been made to improve the uni- 
formity of the lateral deposit of an isolated swath. In actual 
practice, however, uniformity of the deposit, is dependent 
upon not only the distribution of concentration within a 
single swath but also the manner in which the swaths are 
superimposed, that is, the lateral spacing of adjacent passes 
of the airplane. 

Since flight-path spacing is an important factor with re- 
gard to economy of operation, its effect on the total dis- 
tribution of deposit has been determined for several of the 
individual swaths shown in figure 10. These results are 
presented in figure 12. In order to determine the total con- 
centration of deposit, the ordinates of overlapping swaths 
are added as illustrated in figure 12 (a). If an optimum 
deposit is taken to be one in which the minimum required 
concentration is distributed uniformly on the ground, the 
concentration which exceeds this required minimum value 
will constitute an excess of deposit. A ratio of the excess 
deposit (hatched area in fig. 12 (a)) to the total deposit may 
therefore be considered a measure of the uniformity of the 
swath. The variation of this ratio with lateral flight-path 
spacing is plotted in figure 12 (b) for 5„ = 300 microns and 
in figure 12 (c) for 5^=200 microns. The coded lines cor- 
respond to the efflux rates given by figure 10 (a). 

These plots indicate that, for a similar degree of uniformity 
in the swath, the flight-path spacing can be increased by 
modifying the efflux rate in the manner shown in figure 
10 (a). For example, if 10 percent excess deposit is con- 
sidered tolerable, the modified flow rate would permit an 
increase in the flight-path spacing of approximately 0.4 
semispan for 5^=300 microns and 1.0 semispans for 6„=200 
microns. A second conclusion to be drawn from figure 12 
is that, of the two spray spectra considered, the spray with 



Figure 11. — Distribution of deposit with a linear spanwise increase in 
the mean diameter of the droplet-size spectrum. 5„= 200-1- 125?/o: 

C/,= 1.2; Zo = 0.5 semispan; ^=1.0. 

- , Ul/o 


S„=200 microns provides a more uniform swath for equiva- 
lent flight-path spacings. Furthermore, the uniformity of 
this swath appears to be somewhat less sensitive to flight-path 
spacing. This consideration would be of importance in 
cases where it is difficult to estimate the optimum spacing 
between adjacent passes of the airplane. 

CONCLUSIONS 

Results of the foregoing analysis of factors affecting the 
dispersion of spray droplets issued from a moderate-size 
agricultural airplane indicate that: 

1. Greater lateral dispersion of aerial sprays will result 
when the position at which droplets are ejected is moved 
toward the wing tip, when the mean diameter of the droplet 
spectrum is reduced, or when either the airplane lift coeffi- 
cient or the altitude is increased. 

2. The uniformity and effective width of the swath is im- 
proved by increasing the mass effiu.x rate with distance from 
the plane of symmetry. 

3. As compared with a spray spectrum of 300-micron mean 
diameter, the spray spectrum with a mean diameter of 200 
microns allows greater flight-path spacing and produces a 
more uniform deposit. Moreover, the degree of uniformity 
is less sensitive to changes in the spacing of adjacent passes 
of the airplane. 

Langley Aeronautical L.aboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., August 4, 1953. 


<L <t. 

Swath no. I Swath no. 2 



(a) Overlapping of adjacent swaths. 

(b) 5„ = 300 microns. 

(c) 6,„ = 200 microns. 

Figure 12. — Illustration of the excess deposit of overlapping swaths 
and its variation with flight-path spacing. 


APPENDIX A 

EVALUATION OP ERROR IN DROPLET POSITION DUE TO PROXIMITY OF THE BOUND VORTEX 


When a particle is released from the trailing edge of a 
wing, the velocity induced by the bound vortex may repre- 
sent a substantial part of the total velocity acting on the 
particle during the initial part of its motion. In the present 
analysis an attempt to account for this additional down wash 
was made by simply assuming that particles were ejected 
downward from the trailing edge. The initial velocities of 
the particles were assumed to be equal to their steady veloc- 
ities when falling through a uniform stream having an air 
velocity equivalent to the velocity induced by the idealized 
flow pattern at the initial position of the droplet. In this 
appendix, the errors introduced by such an assumption are 
evaluated for the case in which the wing is represented by a 
horseshoe vortex. 

In computing these errors it is convenient to determine an 
air-velocity correction which, when added to the idealized 
flow velocities, will provide the velocities induced by a horse- 
shoe-vorte.x system. The required correction is seen to be 
the velocity induced by a horseshoe vorte.x of span 2b with 
legs extending forward to infinity in the direction of flight. 
The strength of this horseshoe vortex is the same as the 
strength of the vortices in the idealized flow field and its 
sense of rotation is such that, when superimposed upon the 
idealized vortex arrangement, the forward-extending legs 
cancel the vortices which would otherwise e.xtend to infinity 
ahead of the airplane. The coriection velocity induced by 
the horseshoe vortex reflected from the ground plane is 
small in relation to the velocity induced by the bound vortex 
at the wing and in the present appro.ximation will be 
neglected. 

For convenience the origin of coordinates in this case is 
taken at the center of the bound vortex with x positive in a 
rearward direction and z' positive in an upward direction. 
The downwash-velocity correction (denoted by the subscript 
c) in the xy-plane is given by the relation 

{Wa\ _ Cl jir V(l-yy+a:^ . V(l+y)^+? ~|_ 2 | 

\U)c 4vAUL 1—2/ 1+2/ J (1— y)(l+2/)) 

(Al) 

At a:=2.0 semispans, the velocity on the midspan given by 
equation (Al) is approximately 5.8 percent of the downwash 
induced by the idealized flow field when r=0 and the air- 
plane altitude is 0.5 semispan. Consequently, the error in 
droplet position associated with the error in induced down- 
wash beyond a;=2.0 semispans will be assumed to be negli- 
gible. Moreover, it will also be assmned in the following 
calculations that the vortex cores are held fixed during the 
time required for the airplane to travel 2.0 semispans and 
that the downwash-velocity correction does not change with 


the departure of the droplet from the sy-plane. The last 
two assumptions lead to overestimation of the magnitude of 
the downwash correction whereas the first leads to rmder- 
estimation of its duration. 

If the equation for vertical motion of a droplet (eq. (2)) 
in the ideahzed flow field is subtracted from a similar equa- 
tion in which the downwash-velocity correction (eq. (Al)) 
and the increment in drop velocity due to this downwash 
correction have been included, the resulting equation when 
(7 d/?/ 24 is constant is 


\U)c a C'of?r/wA /wA"| (A2) 

dr -R^ 24 LWA W/J 

which may be written 

aCnRdz'c (jCdR (Wa\ 
dr^ '^Ra24: dr ~ Ru24t \Uj, ^ ’ 

An integration of equation (AS) over a distance extending 
from the trailing edge of the wing (x=0.3 semispan) to a 
position 2.0 semispans downstream of the bound vortex gives 
the approximate vertical droplet-position error caused by 
using the idealized flow pattern rather than a horseshoe vortex 
in the neighborhood of the wing. If the horizontal droplet 
position (in semispans) downstream of the bound vortex is 

given by ^-|-0.3, where U is the airplane flight speed, the 

additional downwash (uDalU)^ acting on a droplet released 
from the trailing edge may be determined by substituting 

T—x — 0.3 into equation (Al), where It found 

that over the interval 0.3 ^2.0 semispans and for y values 

equal to or less than 0.75 semispan equation (Al) was essen- 
tially independent of spanwise distance y. An empirical re- 
lation which gives a satisfactory approximation of equation 
(Al) over this interval is 

^ =-4.09-l-6.60T-2.62r' (A4) 


where 0 ^ 1.7. 

The drag-coefficient term 


Cj,R 
24 ' 


which was assumed to be 


constant in equation (A3), wUl now be evaluated for a 
steady-state (terminal droplet velocity) condition. With 


d(—) 

-=~=Q and equation (4) reduces to 

^ "la A (dj)R t Ru bg 

\UJ, “24 ,7 


(A5) 


13 
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where the subscript t denotes a terminal velocity condition 
and this velocity is in a negative direction. Equation (A5) 
may be written in a more convenient form as 


Tt GoRt Ru^ bg 

24 ~ a 


m 


which, together with the functional relation between (7 dS/ 24 
and R (table I), defines CoRtj24, the constant coefficient in 
equation (A3). 

Equation (A3) may now be put into the form 

(A7) 

where 

cr CJl, 

°"~~Ru 24 


vortex, namely, that droplets are ejected at a finite downward 
velocity 



(AlO) 


Let As' be the total difference in distance traveled by a 
particle ejected into an airstream at a finite velocity (eq. 
(AlO)) and a similar particle released from rest into the same 
airstream. If the airstream is considered to have a uniform 
velocity WalU and the drag term OoRtl^i is again taken 
to be constant. As' (in semispans) from the differential 
equation of droplet motion becomes 


As' 



(All) 


The general solution of equation (A7) is 

s'c=Cie”“’'+^(Zo — T^-\-(i2T^-\-G2 (A8) 
With the initial conditions 


r =0 



the arbitrary constants become 

ao Ojl I 2^2 / A 

Cl = — Gi= — 3 (A9) 

a a a 

When equation (A 8 ), with constants defined by equation 
(A9), is evaluated at t= 1.7, the value obtained is the approx- 
imate error in vertical droplet position caused by using the 
idealized flow rather than the flow induced by a horseshoe 
vortex. The downwash is underestimated when the ideal- 
ized flow field is used, so the correction is additive in a down- 
ward direction. 

It remains to determine the influence of the assumption 
used to compensate approximately for neglect of the bound 


The approximate total error in vertical droplet position 
which may be expected as a result of the assumed starting 
conditions used in the trajectory computations is then the 
difference between equations (A8) and (All). This error 
was evaluated at several spanwise starting positions for the 
flight conditions Ci=1.2 and 2o=0.5 semispan and for drop- 
let diameters of 200 and -300 microns (the mean diameters of 
spray spectra to be considered) . 

The calculated errors in vertical droplet position (in 
semispans) for these conditions are as follows: 


ro. 

semisp;ms 

Error* Csemispans) for— 

6=200 microns 

6=300 microns 

0.25 

-0.017 

0. 004 

.50 

-.011 

.013 

.75 

.005 

.053 


*A negative sign indicates that the correction would 
be applied in a downward direction. 


Similar computations were made for the error in lateral 
droplet position caused by neglect of the horseshoe vortex. 
The corrected lateral droplet position was inboard of the cal- 
culated position; however, the magnitude of the correction 
was found to be insignificant. 


APPENDIX B 


LEAST-SQUARES REPRESENTATION OF TRAJECTORY DATA 


In the computation of the distribution of deposit for the 
case of a discrete nozzle the unknown coefficients ki, ^ 2 , and 
ks and the spanwise starting position 2/0 in equation (19) are 
constant. It is desired to determine the values of k^, k^, and 
^3 in the formula 

5=*i+A: 2P+^3P^ (Bl) 

which give the best approximation of n data points (table 
III) at the particular j/o value in question. The principle 
of least squares provides a method of computing these 
coefficients which minimizes the sum of the squares of the 
differences between the values of 5 in table III and the values 
of 5 computed from equation (Bl). Thus the three unknown 


coefficients may be determined by solving the following set 
of simultaneous equations: 


^5=kin-\-k2^P+k3^P^ 

'y',P5=k,y^,p+k2'>^,p^+k^y^,p^ >- 

^P^5=k,-^P^-kk2^P^+k,J2P^J 


(B 2 ) 


The deposit for discrete nozzles was evaluated only at the 
yo positions for which trajectory computations were available. 
This limitation obviated the necessity of interpolating to 
determine the coefficients at other positions; however, in 
determining the deposit for the case of a continuous spanwise 
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distribution of mass flow (eq. (16)), the variations of ki, k 2 , 
and kz with must be known. From experience gained 
while determining the coefiicients of equation (B2) for 
various values of y^, it was found that the variation of these 
coefficients with yo could be satisfactorily approximated by 
the relations 

*i=a22/o^+a3l/on 


*2=&li/o + &2j/o^ 


(B3) 


fc3 = Cli/0 + C2?/o^ J 


6i=216.9 
62 = -49.76 
c,= -0.0537 
C2=-8.141 


The comparison of the points determined by trajectory 
computations and the empirical representation of these 
points by equation (B4) is shown in figure 13. 


Equation (Bl) then becomes 

5 = a^ya^ + fflgi/o® + hiyoP + hiy^^P + + c^y^^P^ (B4) 

Equation (B4) contains six unknown coefficients, and for the 
case in which 6i=1.2 and 2 q= 0.5 semispan (the flight 
condition for which the deposit with a continuous mass 
flow was evaluated) there is a total of 20 data points (see 
table III). Applying the method of least squares to equa- 
tion (B4) provides a set of six simultaneous equations: 

= o,2^ya + + b 1 ~^^Pyo + b^ yo"* + 

ci12PV +c2j2P''yo* 

Y^yo^5= a2Sf/o" + O 3 + biYlPya + bi ^Pyo^ + 

Cl + csZF^yo® 

Z-Pyo5 = tti'^Pya^ + a^^Pya^ -1- b,J2P%^ + b2^P^ -|- 

c 1 Z-P^o'^ + Ci'^Phjo^ 

J^Pyo^d = a^ZFyo^ + asZFyo^ + b.J^P^y,^ + b^lZP^y ,^ + 

ciZFV+c2ZP“V 

ZF^y„5 = a,JZPV + a,-^Phj,^ + h ZF^2/o^+ &2ZF^l/o^+ 

Ci Z^Vo^ + C 2 ZF^o^ 

ZF^yo^S = U2l2P%* + a,J2P^yo^ + b,J2P^y,^ + f.2ZF^2/o^+ 
ciZFV+c2ZFVo'* 

These equations may be solved for the unknown coefiicients 
(ffl 2 , « 3 ) and so forth) by any convenient method of solving 
linear simultaneous algebraic equations. 

The following numerical values of the coefficients in 
equation (B4) were determined for the case in which 6/,= 1.2 
and 2 q= 0.5 semispan: 

®2=373.7 

«3=- 157.4 



Fiourb 13. — Empirical representation of droplet diameter as a func- 
tion of yo and y^. C;,= 1.2; Zo — 0.5. 
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TABLE I— VALUES OF C!dE/ 24 AS A FUNCTION OF JJf 


i? 

Cz)i?/24 

R 

CdRI^A 

0 

1.00 

14.0 

2.009 

.05 

1.009 

16.0 

2. 109 

.1 

1.018 

18.0 

2. 198 

.2 

1. 037 

20.0 

2. 291 

.4 

1. 073 

25.0 

2. 489 

.6 

1. 103 

30.0 

2. 673 

.8 

1.142 

35.0 

2. 851 

1.0 

1. 176 

40.0 

3.013 

1.2 

1.201 

50.0 

3.327 

1.4 

1.225 

60.0 

3.60 

1.6 

1.248 

80.0 

4.11 

1.8 

1. 267 

100.0 

4. 59 

2.0 

1.285 

120.0 

5.01 

2.5 

1.332 

140.0 

5. 40 

3.0 

1. 374 

160.0 

5.76 

3.5 

1.412 

180.0 

6. 16 

4.0 

1.447 

200.0 

6. 52 

5.0 

1.513 

250.0 

7. 38 

6.0 

1.672 

300.0 

8. 26 

8.0 

1.678 

350.0 

9.00 

10.0 

1.782 

400.0 

9. 82 

12.0 

1.901 

500.0 

11.46 


tFroni reference 3. 


TABLE II. — AIRPLANE AND SPRAY CHARACTERISTICS 


Wing semispan, 6, ft 

Wing aspect ratio, A. 

Density of droplets, pd, slugs/cu ft 

Density of air, pa, slugs/cu ft 

Lift coefficient of airplane, Cl 

Flight velocity, C , ft/sec 


19.5 
5. 35 
1 55 
0. 002375 
1.2; 2.0 
85. 0; 62. 6 


TABLE III.— RESULTS OF DROPLET-TRAJECTORY 
COMPUTATIONS 


yo, semi- 
spans 

Vg, semi- 
spans (*) 

8, mi- 
crons 

Vdg, ft/sec 

Wd,,itlsec 

tg, sec 


(a) Cr,=1.2; zo 

=0.5 semispan 


0. 25 

0. 506 

210 

2. 09 

-1.97 

2. 88 

.25 

.441 

275 

2. 08 

-3. 08 

2. 26 

.25 

.385 

375 

2. 06 

-4. 21 

1.71 

.25 

.347 

500 

2.01 

-5. 65 

1.31 

.25 

.313 

700 

1. 83 

-8. 75 

.95 

.50 

1.675 

150 

10. 38 

-1. 11 

3. 03 

.50 

1.130 

210 

9. 95 

-2. 12 

2. 18 

.50 

.925 

275 

7. 84 

-3. 03 

1.77 

.50 

.795 

375 

6. 45 

-4. 40 

1.42 

.50 

.710 

500 

5. 13 

-6.83 

1. 15 

.50 

.628 

700 

4. 66 

-9. 16 

.84 







.625 

1.670 

210 

8. 76 

-i.67 

2. 42 

.625 

1.250 

275 

11. 47 

-2.91 

1.61 

.625 

1.006 

375 

10. 03 

-4. 44 

1.35 

.625 

.893 

500 

8. 93 

-5.71 

1.00 

.625 

.787 

700 

6.77 

-9. 51 

. 75 







.75 

1.77 

275 

6. 56 

-2. 56 

2. 26 

.75 

1.33 

375 

11.03 

-3.84 

1. 25 

.75 

1. 10 

500 

11.01 

-6. 02 

. 95 

.75 

.95 

700 

8.82 

-9.66 

.65 


0>) 

Ci~2.2; 20=0.5 semispan 


0.25 

0. 692 

150 

2. 56 

-1.36 

3.72 

.25 

.556 

210 

2.71 

-2. 06 

2.61 

.25 

.480 

275 

2. 76 

-2. 98 

2. 03 

.25 

.415 

375 

2. 76 

-4. 26 

1.54 

.25 

.370 

500 

2. 72 

-5.97 

1. 23 

.25 

.330 

700 

2. 46 

-8.91 

.94 

.50 

(••) 

150 




.50 

1.290 

210 

i5. 14 

-2. 45 

1.86 

.50 

1.025 

275 

11.51 

-3. 16 

1.52 

.50 

.855 

375 

8.79 

-4. 53 

1.25 

.50 

.750 

500 

7.64 

-5. 67 

1.02 

. 50 

.653 

700 

6. 26 

-9. 64 

.76 

.625 

2. 305 

210 

6. 16 

-1.93 

3.51 

.625 

1.430 

275 

15. 00 

-2. 68 

1.45 

.625 

1.095 

375 

13. 90 

-4. 54 

1.06 

.625 

.940 

500 

11.60 

-6. 57 

.87 

.625 

.815 

700 

9.11 

-10. 10 

.67 







.75 

2.210 

275 

5. 94 

-2. 42 

3.11 

.75 

1.475 

375 

14. 03 

-3. 69 

1. 18 

.75 

1. 155 

500 

14. 92 

-6. 23 

.75 

.75 

.985 

700 

12. 64 

-10. 40 

.65 


(0 

Cx = 1.2; zo=1.0 semispan 


0. 25 

1.750 

150 

6. 56 

-1. 13 

6. 55 

.25 

.851 

210 

5. 36 

-2. 12 

4. 47 

.25 

.500 

375 

3.01 

-4. 21 

2. 35 

.25 

.360 

700 

1.97 

-8.81 

1.75 







.50 

1.925 

210 

4. 16 

-1.90 

5. 60 

.50 

1.070 

375 

6. 29 

-4. 08 

2. 74 

.50 

. 740 

700 

4. 21 

-8. 58 

1.70 


♦The Ve values are without cross wind. To account for a cross wind, the product of tg and 
the cross-wind velocity component may be applied as a correction to the tabulated i/x values. 
♦•Computations were discontinued before droplet intersected ground plane. 
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